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In this paper we consider networks of dynamical systems that evolve in synchrony and investigate 
how dynamical information from the synchronization dynamics can be effectively used to learn the 
network topology, i.e., identify the time evolution of the couplings between the network nodes. To 
this aim, we present an adaptive strategy that, based on a potential that the network systems seek 
to minimize in order to maintain synchronization, can be successfully applied to identify the time 
evolution of the network from limited information. This strategy takes advantage of the properties 
of synchronism of chaos and of the presence of different communication delays over the network 
links. As a motivating example we consider a network of sensors surveying an area, in which 
information regarding the time evolution of the network connections can be used, e.g., to detect 
changes taking place within the area. We propose two different setups for our strategy. In the first 
one, synchronization has to be achieved at each node (as well as the identification of the couplings 
over the network links), based solely on a single scalar signal representing a superposition of signals 
from the other nodes in the network. In the second one, we incorporate an additional node, termed 
the maestro, having the function of maintaining network synchronization. We will see that when 
such an arrangement is realized, it will become possible to effectively identify the time evolution of 
networks that are much larger than would be possible in the absence of a maestro. 



I. INTRODUCTION 



In recent years, much effort has been devoted to the study of synchronization of networks of coupled dynamical 



systems 
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been obtained under the assumption that communication between the connected systems is immediate; i.e., the signal 
from system i at time t is received by systems j at the same time t. Most of these studies have also focused on the 
case where the network topology is static, i.e., the coupling strengths between the connected systems is fixed and 
does not evolve in time. Yet in real situations networks of dynamical systems are often time-varying and affected 



by communication delays. In this paper we will address situations incorporating both the time varying structure 
common in typical real networks and the delays which inherently affect communication across the network. Our main 
emphasis will be on the situation where the time-dependent connection strength from system j to another system i 
in the network is unknown by system i, and we investigate how system i can deduce these strengths by an adaptive 
strategy making use of the phenomenon of synchronization of chaos. As we shall show, diversity in the delay time 
associated with different links, as well as chaos of the exchanged signals, will play a crucial role in achieving this goal. 

As a possible motivation, we will consider as our reference application that of a network of sensors that are dispersed 
over a ground and interact with each other through direct wireless communication. Namely, we will assume that the 
sensor network is employed to detect the evolution of the configuration within the surveyed area, e.g., due to motions 
of objects on the ground. In that case, a movement of an object across a communication pathway affects the strength 
of the signals exchanged by the nodes at its endpoints. Therefore, if we can succeed in identifying the time-evolution 
of the network couplings, this will provide information potentially enabling us to localize the position of the objects 
moving across the ground. We consider that each of the dynamical systems ('nodes') are identical chaotic systems 
broadcasting chaotic signals. As the network evolves, each node uses an adaptive strategy (which we specify) to 
maintain synchronization, and, in so doing, it determines the strengths of the couplings to it. 

Moreover, we will also show how our adaptive strategy could greatly benefit from the introduction of a special node, 
having the specific function of maintaining network synchronization. In what follows, we will refer to this particular 
node as the maestro, in order to distinguish it from the remaining nodes of the network, which we will refer to as 
the orchestra. For example, for the case of our reference application, we will make the assumption that this node is 
positioned in such a way that the strength of the signal sent by this node is not affected by the movements of the 
objects on the ground (e.g., it might be located atop a tower, hill or tall building). We will see that when such an 
arrangement is realized, it will become possible to effectively identify the time evolution of much larger networks than 
would be possible in the absence of a maestro. 

The problem of identifying the structure of real complex networks, arises in all those situations where the exact 
network topological structure is unknown or uncertain, therefore it is of relevance in various fields of science and 
technologies. Examples are metabolic networks, biological neural networks, electric power grids, and so on 
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Moreover, in many real world situations, there is often an intrinsic difficulty in extracting information on the network 
structure. Hence, the problem of learning the network top ology, re presents a challenge in a number of real applications, 



that involve both natural and man-made systems 
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In Sec. II we present our proposed adaptive strategy for maintaining synchronization of chaos and determining 
network topology. We do this both for the case with no maestro and for the case with a maestro. We also show that 
it is essential for our technique that signals be chaotic and experience time delays along the links of the network. In 
Sec. Ill we test our strategy through numerical experiments. We show that the presence of a maestro can greatly 
benefit our strategy. Conclusions and further discussion are presented in Sec. IV. 



II. ADAPTIVE STRATEGY 



In a previous paper 29( we showed how a time-varying complex network can be synchronized by using a 'potential' 
that each network node seeks to minimize. In this paper, we will extend this concept with the goal, of not only 
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maintaining synchronization, but of also identifying the time-varying evolution of the network itself. In Refs. 
it was shown how to identify the static structure of a complex network of synchronizing systems by introducing a 
response (replica) network that can adjust its couplings to converge to those of the original one. The main focus of our 
paper will be on identifying the evolution of time-varying couplings of a complex network from limited information. 
We will show how this may be achieved by solving a small number of auxiliary differential equation, provided that 
some general conditions on the time delays of the network links are met. We will consider two versions of our problem. 
In Version I the maestro is absent, while in Version II we have a maestro. 

In Version I, the set of dynamical equations describing the dynamical state vector Xi(t) of dimension n at each 
node i of the time evolving network is the following: 

N 

x i {t)=F(x i (t))+T{r i {t)- J2 Aj^Hixiit-Tij))}, z = l,2,...,V, (1) 

where 

N 

n(t)= £ A l3 {t)H{x 3 {t- nj )) (2) 

3=1 ,j¥=i 

is the incoming signal at node i, which is a linear combination of the different scalar transmitted signals H(xj(t)), 
through the time varying coupling coefficients Aij{t). Here T = [T\, I?2, r„] T is a constant n-vector specifying 
the coupling of the scalar received signal to node i. (In our numerical experiments in Sec. Ill, we will take the 
components of T to be all zero except for one component, say £, for which Yg > 0.) We assume that the time delays 
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Tij, between each pair of connected nodes i and j, are known, e.g., the sensors occupy fixed positions and thus the 
communication delays can be easily obtained from knowledge of the relative distances between the sensors. In contrast, 
we assume that the strengths of the network interactions, represented by the Ay, are unknown and variable in time. 
Here, the Ay(t)'s might represent the spreading, scattering, and attenuation of a wireless output signal broadcast by 
node j to node i. In the case where the input to i from j is accomplished by propagation of an electromagnetic wave 
from j to i, the link strength Ay will vary with time as a body moves through the propagation path from j to i. 
Although the Ay(i) vary with time, we shall be most interested in the case where the time scale for this variation is 
slow compared to the natural time scale of the chaotic dynamics of an individual uncoupled system, 

±i(t) =F( Xi (t)). (3) 

The quantity Ay , appearing in Eq. ([1]) represents an estimate made at node i of the true coupling Ay (t) . How this 
estimate can be made, as well as the accuracy of the estimate, will be the main focus of this paper. 

In Version II, we consider an additional dynamical state vector, i.e., XM(t), associated with the maestro node, and 
the following set of dynamical equations: 

N 

±i(t) = F(xi(t))+T{ n -[j^j A l0 H{x t {t-T lj )) + {l-a)A l M{t)H{x l {t~T lM ))]}, i = 1, 2, N, (4) 

where 

N 

n ^ = N~T ^ A ij {t)H{x j {t-T ij )) + {l-a)A i M{t)H{x M {t-T i M)), (5) 

and the maestro obeys the same dynamical evolution as an uncoupled node, i.e., 

x M (t) =F(x M {t)). (6) 

Here AiM(t) is the strength of the coupling of the maestro to each node i = 1, N; tim is the communication delay 
between the maestro and node i = 1, N; 0<a< 1 is a parameter that tunes the relative strength of the coupling 
exerted on the network nodes by the maestro, with respect to that of the other nodes in the network (the orchestra). 
When a = 0, each node is affected only by the forcing by the maestro, and therefore it can be seen as an isolated 



slave system connected to the master. On the other hand, when a = 1, each node will not feel the influence of the 
maestro and synchronization can only be realized from the mutual coupling with the other nodes in the network (the 
orchestra); thus the case a — 1 corresponds to our problem in Version I. Values of a between and 1 (the case of 
most interest) correspond to intermediate situations where each node feels both the influence of the maestro, as well 
as that of the orchestra. 

In what follows we assume that Aim if) is known a priori at each node i. In Sec. III.D we also address the case 
that the assumed known coupling with the maestro is affected by a small error. 

In the particular case where we actually have the correct estimates Aij — Aij, the network equations in ([I]) and 
(H|), can be rewritten as 

N 

±i(t) =F(xi{t))+T H(x i (t-T ij ))}, i = l,2 t ...,N, for Version 1, (7a) 

and 

N 



± i (t)=F(x i {t))+r{j^— [ Yl A lJ {t){H(x J {t-T l] ))-H(x l {t-T lJ ))}+ 



a 

T 

+(l-a)A i M(t){H(x M (t-T i M))-H(x i (t-Tm))}}, i = l,2,...,N, for Version 2. (7b) 



Thus, in a similar way to the case of a network of coupled oscillators without delays or time variability of the coupling 
Aij a], the synchronization manifold, 

Xi = x-2 = ... = ccjv, for Version 1, (8a) 
x\ = X2 = ... = xn = ijf, for Version 2, (8b) 



represents an invariant set for the set of equations (f7aj) or (I7bl) . 

Once the dynamical functions F and H are assigned, the stability of synchronized states to perturbations from 
the synchronization manifold will depend on the dynamical evolution of the network adjacency matrix A(t), on the 
overall coupling T, and on the time delays Tjj (plus on AiMif) and t^m in the case of Version 2). In what follows, we 
will proceed under the assumption that the synchronization manifold (|8al - I8b|) is stable for our choices of the time 
evolution of A(t), Ai^(t), T and the time delays r»j and t,m, and we will focus on the problem of identifying the 
network structure (i.e., obtaining the estimates Ay(i)). 



We assume that at each node i direct information on the couplings A y - (t) is unavailable. Thus, the estimates Ajj at 
node i must be obtained solely using the aggregate incoming signal fj(t). In order to make this estimate, we introduce 
a 'potential' 'I'j at each node i: 



*i(*) = * / e- 1 '( i - t '){r i (t / ) - -j^Y. Au(t)H{xi(!? - r tf )) - (1 - a)A iM (t)H( Xi (t' - r iM )} 2 dt', 



for Version 1, 
(9a) 

for Version 2. 
(9b) 



If the time scale for variation of Aij is larger than v , then ^i(t) can be approximated by 

ft 

*i(t) e-^ t -*'){r l (t') - £ AijWHixiit' - r y ))} 2 < for Version 1, 

(10a) 

/t 
e -^(i-t'){ n (^) _ ^y^Z Aii(t')H(xi{t' - Tij )) - (1 - cOAmWO^' - W)} 2 *', for Version 2. 

(10b) 



As given by (jlOap and (|10b[) . the potential ^i is approximately the time averaged squared synchronization er- 
ror, and v^ 1 is the temporal extent over which the averaging is performed. Note from (jlOap that \&i(t) = 
if £#i^j - T i:j )) = V ; ,, - r (i )). Similarly, from (lObJ we note that $ 4 (t) = Q if 

J2j 9 aA-i 3 (t)H(x i (t - T^)) = J2jjti-Aij(t)H(xj(t - nj)) and - t iM ) = - t,m)- Due to the chaotic na- 

ture of the x's, provided that the Tij are all distinct, these conditions can only be satisfied if Ay(f) = Ajj(t) (i.e., 
our estimates of the time evolutions of the network couplings are correct) and Xi(t — r^) = Xj(t — Tij) (i.e., the 
network systems are synchronized in time). Thus we seek to minimize ^i. Therefore, we impose that the Ajj(t) 
evolve according to the following gradient descent relation: 



" "/3-TT 1 = ~ X! AkCijk), for Version 1, 



dt 



dA 



d AlL - -Q^L - 2-^-\c- ■ + (1 - a)£ 



N ■ 



-j ^ AjkCjjk], for Version 2 , 



(11a) 
(lib) 
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where 



Cij(t) = v / e-"t*-*'V i (t , )fl-(a! < (t / - nj))dt\ (12) 



n 



C,ifc(*) = v ! e-^-^Hixiif - nj))H(xi(tf - n k ))dt', (13) 



o 



&i(t) = f / e-'C-'O^^ff^lf _ Tij^Hixii? - r iM ))dt', (14) 



and /3 is a constant larger than zero. 

In practice it is useful to avoid doing the above integrations at each time step of our adaptive computation. Rather, 
we note that Cjj(i), Cijk(t), and £jj(f) given by (fT2"|) and (fT^j) . and (TT4")) satisfy simple differential equations, and we 
employ these in our adaptive computation of Aij(t): 



dcjj(t) 



vdjit) + vn(t)H(xi(t - Ty)), (15) 



dC 

= -vC ijk {t) + uH( Xi (t - nj^Hfait - r lk )), (16) 



^ = -v&j(t) + vA m (t)H{xi(t - nj ))H{xi{t - t iM ))- (17) 



The simplest version of this setup, used in our numerical example in Sec. Ill, corresponds to letting (3 — ► oo, in 
which case Eqs. (11) implies that we can obtain the Aik by solving the system of linear equations, 



Cij = ^^AikCijk, for Version 1, (18a) 

k=£i 

cy - (1 - a)iij = — _ - 2^ A ikCij k for Version 2, (18b) 
where cy, Cijk, and £jj are obtained by solving (fl5|) . (fTBI) . and (fT7|) . 
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According to our previous discussion (in particular, the condition for the approximate equality in Eq. (I10a|) ^l . in 
order for our identification strategy to work properly we should choose v such that it is in the intermediate range, 

T s < v- 1 < T n , (19) 

where T s is the time scale of the chaotic dynamics of an uncoupled nodal system (x — F(x)), and T n is the time scale 
on which the network itself changes (i.e. the time scale on which Aij(t) varies). We will subsequently give quantitative 
definitions of T s and T n for our examples. 

We emphasize that we try to identify the time evolutions of the several Aij's, by using only one piece of external 
information at each node i, namely, the incoming signal r^t). Our ability to extract estimates for the N — 1 values 
of Aij(t) (i — 1,2, ...,N;j ^ i) at each node i from each incoming signal, relies on the temporal signal diversity of 
Ti(t) due to the presence of both chaos and different time delays along the network links. To illustrate the necessity 
of different time delays, consider the following example. Assume, for instance, that our network consists of three 
coupled systems in the Version I scenario (N = 3 in Eq. (1)), and we want to identify the time evolutions of Ai2(t) 
and A\3 (t) when the only available information at node 1 is r\ (t) , the time evolution of the incoming signal at node 
1. If we are on the synchronization manifold, then Eqs. (1) and (2) yield 

x x {t) = F( Xl (t)) + T{[A 12 (t) - A 12 (t)]H( Xl (t - r 12 )) + [A 13 (t) - A 13 (t)]H( Xl (t - r 13 ))}, (20) 

and, for Tyi = T13, the H terms are equal, yielding for the term in (f2T))) that is in curly brackets 

H( Xl (t - raa)){[Ai 2 (t) - A 12 (t)] - [A 13 (t) - A 13 (t)}}. (21) 

Thus, if T12 = T13, we see that we could obtain synchronization at node 1 for any A\2(t) and Ai 3 (t) such that 
Ai2(t) + Ai 3 (t) = Ai2(t) + Ai 3 (t), and we cannot hope to estimate both Ai2(t) and Ai 3 (t) only from ri(t). In 
contrast, if T12 ^ T13, then on the synchronization manifold x 2 {t — T12) =/= x 3 (t — T13) and the previous degeneracy is 
removed. Thus we now have the possibility that it may be feasible to obtain the desired two unknowns Ai2(t) and 
Ai 3 (t); in other words, it may become possible to obtain good estimates Ai2(t) and A\ 3 (t), provided that t\2 7^ Tn- 
Hence the presence of different delays is crucial in the identification process. 

Moreover our strategy benefits from the chaotic nature of the exchanged signals at the network nodes. In particular, 



9 

we take advantage of the high information content encoded in chaotic signals. As an illustrative counterexample, 
assume that, in place of the chaotic signals, we had sine waves of the same frequency broadcast by the network 
systems. The communication delays would result in phase shifts of the sine waves, and the signal received at each 
node, being the sum of sinusoids of the same frequency, would be a sinusoidal signal. In this case, there would be 
no information content in the received signal other than its amplitude and phase; therefore, we would not be able to 
extract more than two pieces of information from it. 



III. NUMERICAL EXPERIMENTS 



A. Version I 



As a first example, we consider the problem in Version I [Eqs. (JlJ and |J2J], and we seek to estimate the time 
evolution of a small network of N — 3 nodes and M = 6 directed links (two directed links incoming to each node) . 
For our uncoupled dynamical system, ±i = F(xi), we consider a Rossler oscillator, xt = {xn,Xi2,x a ) T and 



F( Xi ) 



-x i2 - x a 
xn + 0.165x i2 
0.2+ (2*1 - I0)x a 



We choose the oscillators to be coupled in the xn variable, i.e., H{x) — x\, V — [ri,0, 0] T , Ti > 0. 
The delays associated with links j — > i, are given in matrix form, T = {t^}, as follows: 



(22) 



T = 0.3T, 



1 2 

1 3 

2 3 



(23) 



where T s is defined as the time s at which the autocorrelation function C{s) of x s \{t), C(s) =< x s i(t + s)x s i(t) > — < 
x s i(t) > 2 , decays to C(0)/2, and x s (t) is a typical orbit on the chaotic attractor of the uncoupled system, x s = F(x s ). 
Note that for each i, the r^- for every j ^ i are all unequal. For our choice in Eq. (|3H)) . we obtain T s = 0.75. 
For t > we assume the following network evolution, 



Aij(t) = 1 + dj sin(u>ijt + 



(24) 
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where the e.^ = tji are random numbers drawn from a uniform distribution between and 1, the u>ij = ojji are random 
numbers drawn from a uniform distribution between u> m in and Lo max , and the = (jP^ are random numbers drawn 
from a uniform distribution between and 2tt. We define the network evolution time scale T n for our example to be 

1 n — 10 max • 

The network state vectors are initialized as follows, 



x n = x i + cpie ix , x° l2 = x 2 + cp 2 e ly , x% = x° + cp 3 \e iz \. (25) 

Here (x\, x%, x®) is a randomly chosen point on the Rossler attractor, obtained by evolving the uncoupled system, 
x s = F(x s ), until it is on the attractor, and then taking its state at some arbitrarily chosen subsequent random 
time; e^, ti V and et z are zero- mean independent random numbers of unit variance drawn from a normal distribution; 
pi = 7.45, p2 — 7.08, P3 = 4.25 are the standard deviations of the time evolutions of x s i(t), x S 2(t), x S 3(t) obtained 
from numerical solution of an uncoupled oscillator, x s — F(x s ). For our numerical experiment we take Ti = 0.25, 
v = ~ 0.17, w max = 0.4 x 10~ 4 and Lo m in = w maa; /2. Using these parameters we have that the inequalities in Eq. 
(UHJ) are satisfied: 0.75 < 6 < 2.5 x 10 4 . The computations were done over a time interval < t < 10 5 . 
To measure the extent to which our strategy works we introduce the following two error measures: 



E Xl (t) = -J— V \xa(t)- < xn(t) > |, (26) 

where < xn(t) >= ^St^iiW 

Figures [TT21 and[3]show results of our computations. Figure Q] shows xn(t) (Fig. [T{a)) and Xii(t) (Fig. [TJb)) versus 
t at the end of the the run, 0.999 x 10 5 < t < 10 5 . For all nodes, i = 1,2, 3, xu(t) and X2i(t) are identical to within 
the width of the plotted curves. Figure HJa) shows that EA(t) is typically less than 0.5% over the entire run, while 
Fig. [Hb) shows that the synchronization error E Xl (t) is typically less than 0.2%. These results confirm that our 
adaptive strategy is effective in synchronizing the network and that the j4y (t) closely follow the evolutions of the time 
evolution of the true couplings Aij (t) (see Fig. \S§ . 
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FIG. 1: Time evolutions of Xu(t),Xi2(t) between t = 99,900 and t = 100,000. The plots for i = 1,2,3 are identical to within 
the width of the plotted lines. 



FIG. 2: Time evolutions of EA(i) and E xi (t), for t between and 10 5 



Figures [4] [5l and [6] show the results of another simulation where we considered a network of N = 5 nodes and 



4 5 1 2 
4 12 3 

T = (UT'. 5 1 3 4, (*si 

1 2 3 5 

2 3 4 5 

with all the other conditions the same as those used in obtaining Figs. 1-3. As can be seen from Fig. [5j for this 
case our identification strategy fails, with the deviating from the evolution of the Aij. At the same time, we 
observe from Fig. [4] that the network systems evolve approximately synchronously in time. Thus, in this case we can 
approximately synchronize the network, but we fail in correctly estimating the evolution of the network couplings. 

Unfortunately, by using the strategy in Version I, we were not able to synchronize complete graphs of Rossler 
oscillators, as the network dimension N was increased. As we will explain in the Appendix of this paper, we have 
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FIG. 3: (Color online) The time evolutions of Aij(t) and Aij(t) are shown from t = to t = 10 5 and agree to within the widths 
of the plotted lines. Each subplot represents a directed link in the network: (a) i=l, j=2; (b) i=2, j=l; (c) i=3; j=l; (d) i=l, 
j=3; (e) i=2, j=3; (f) i=3 j=2. v = §i- ~ 0.17, w m ai = 0.4 x 10" 4 and Umin = Wmm/2. 
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FIG. 4: Time evolutions of a!ii(t),a;«(*) between t = 99,900 and t = 100,000. The plots for i = 1,2,3,4,5 show almost- 
synchronized behavior. 



an understanding of the reason for this failure. Briefly, a requirement for correctly estimating Ay at any given node, 
say i, is that the ry's associated with i are different enough (we address this issue more rigorously in the Appendix). 
However, if N is too large, we cannot make the Ty too separated, because stability of the synchronization manifold 
(ISajl (under the condition Ay = Ay), is lost if the ry's are too large. Thus, for our setup in Version I, achieving the 
required diversity in the Ty's, contrasts with our requirement of having the network synchronize. 

We have also tested our Version I scenario for a case as in Figs. 1-3, but with the network evolving on a faster 
time scale Lo max = 2 x 10~ 4 (instead of 0.4 x 10 -4 ), and for this faster evolution we have found that our strategy fails 
(results not shown). 
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FIG. 5: Time evolutions of EA{t) and E xi (t), for t between and 10 5 
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FIG. 6: (Color) Each subplot represents a directed link in the network. The time evolutions of Aij{t) (black) and Aij(t) (red) 
are plotted from t = to t = 10 5 , with v — gT- ~ 0.17, ui max — 0.4 x 10~ 4 and ui min — uj m ax/2- 



B. Version II 



The failure of our method with large N when applied to Version I is the main reason for considering in what follows 
the alternative set up of our problem in Version II [Eqs. ((4]) and (|5|)], which incorporates a maestro node, having the 
function of maintaining the network synchronizability. Thus we assume that each node/dynamical system is subject 
to two actions: a coupling exerted by the other nodes of the network, which in what follows we will refer to as the 
orchestra and a coupling exerted by a master system, which in what follows we will refer to as the maestro. By tuning 
the parameter a in (j4|) and ([5|) between and 1 we are able to interpolate between a situation where the nodes feel 
only the influence of the maestro and a situation where they feel only the influence of the orchestra (corresponding 
to Version I). 
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We assume that the maestro coupling Aim is known and time independent, but that the other network couplings, 
Aij, are time- varying and unknown at the network nodes. In order to guarantee that the synchronization manifold is 
stable, we take the time delays with the maestro t%m to be small enough. Our hope is that, once synchronization is 
ensured through the maestro, we will potentially be able to choose the other Ty's so that they are sufficiently spread 
that we have a much greater likelihood of being able to correctly estimate the Ay(i). As we will show, by making use 
of the maestro node, we will be able to identify larger networks, whose time-evolution is faster than would be possible 
in Version I (u> max > 10 ). 

To test our strategy in Version II, Eq. (j4|), again, as for Version I, we assume that we want to identify the time 
evolution of globally connected networks of Rossler oscillators, with H(x) = x%, T = \Tx, 0, 0], T% > 0. For our first 
experiment we take a network of N = 5 nodes and M = 20 directed links to be identified. We assume t^m =0.1, 
AiM(t) = 1 at each nodes i in ((H); we have verified beforehand that in the case where Ay = Ay for i,j = 1,...,N, for 
this choice of the parameters, the synchronization manifold (|8b[) is stable. We are interested in identifying the time 
evolutions of all the Ay's associated with the network links. 

We take the time delays Ty along the network directed links i — ► j to be 

nj = T iM + [1 + mod(i + j, N)]9T S , % ± j (29) 

where mod(a, b) is the remainder of the integer division of a by &, and 9 is a variable parameter to be specified. Note 
that, for our choice of the time delays in (|2"9"|) , we have that ry = Tji, i ^ j. Moreover, from (f2"9"| we also have that 
for any node i and for any pair of delays associated with any two links incoming to i, say and r^, j ^ k we have 
that 7^ t^, and the values and Tik are separated by an integer multiple of 0T S . Thus, we do not have equal 
delays along incoming links to any node. 

A first experiment involving our strategy in Version II, is shown in Figures [7| and [HJ For this experiment, we took 
a = 0.2 (i.e., synchronization depends mostly on the coupling with the maestro), and set 9 — 0.6, indicating that 
the separation time between any pair of network delays is a multiple of 0.45 = 0.6T S . Moreover we assumed that the 
network evolves much faster than in the cases reported in Figs. 1-6 of Sec. III. A for Version I, that is, io max = 4 x 10 -3 , 
and uJmin = uj max /2. The network was evolved for a long time until time 2 x 10 4 , starting from synchronous initial 
conditions and assuming Aij(0) = Aj_y(0). 

As shown in Fig. [71 EA(t) < 0.05 and E Xl (i) < 4 x 10~ 3 , throughout the length of the run. The black dots in Fig. 
[8] show the time evolutions of the Ay's and the red lines those of the Ay's. As can be seen, our adaptive strategy 
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is effective in synchronizing the network with the Aij(t) closely following the evolutions of the Aij(t). Remarkably, 
we have succeeded in identifying the time evolutions of all the M — 20 links present in the network. The same 
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FIG. 7: Time evolutions of E A (t) and E X1 (t), for t between t = and t = 2 x 10 4 , a = 0.2, v = 8T S _1 ~ 0.17, 9 = 0.6, Ti = 1, 
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FIG. 8: (Color) Each subplot represents a directed link in the network. The time evolutions of Aij(t) (in black dots) and Aij(t) 
(in red curves) are plotted from f = to f = 2 x 10 4 . 



experiment in Figs. [7] and [8] has also been repeated for different values of a. In Fig. [9j the time averaged errors 
< E X1 > and < Ea > are plotted as a is increased between 0.005 and 0.95 (continuous line). As can be seen, our 
strategy is observed to be effective in identifying the network evolution for values of a up to about 0.45. 

We tested our strategy for increasing network size N up to TV = 18. As shown in Fig. [TQ1 we found that our 
strategy was effective in identifying the network evolution for 5 < N < 18. Note that for our choice of fully connected 
networks, the total number of unknowns to be identified is equal to N(N — 1), e.g., it is equal to 20 when N = 5 
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FIG. 9: (Continuous line) The time averaged errors < E X1 > and < Ea > are plotted versus a, for N = 5, Ti = 1, = 0.6, 
v = 1/(8T 3 ), tOmax = 4 x 10~ 4 , L0 m i n = io max /2. (Dashed line) The same experiment has been repeated for a case where a 
small error in the coupling with the maestro is present, i.e., H = 0.02 (see Sec. III.D). 

and is equal to 306 when N — 18. We expect that, the task of identifying the evolution of the couplings becomes 
more difficult as the network dimension TV is increased. Consider for example a globally connected network of N = 18 
nodes. In such a case, each one of the nodes has to identify the time evolutions of (N — 1) = 17 different couplings to 
the remaining nodes in the network, from only one available received signal. Note that, by making the assumption of 
dealing with fully connected networks, and for which all the links are time-varying, we are considering a worst case 
scenario; in the case of our sensor application (see Sec. I), we expect that our strategy could benefit from a situation 
in which only a fraction of the network links might be changing as an object crosses the surveyed region. 

We also tried to estimate the time evolution of globally connected networks whose dimension N is larger than 18. 
For an TV = 20 node network, our strategy failed at identifying the couplings evolution; however, we repeated the 
same experiment with N = 20, but with the network evolving more slowly (i.e., uj max = 10~ 4 ), and in this case we 
were able to correctly identify the time evolving couplings. 




FIG. 10: The time averaged errors < E X1 > and < Ea > versus the network size JV for 5 < N < 18, with a = 0.2, Ti = 1, 
8 — 0.6, v — 1/(8T S ), uj m ax = 2 x 10~ 4 , LOmin = The crosses indicate the values of < E X1 > and < Ea > obtained for 

N = 5, 6, 7, 8, 9, 10, 12, 15, and, 18. 
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C. Sensitivity to non-identicality of the network systems 

In practical situations (including our reference application of sensor networks) it is impossible to guarantee that all 
the network systems are precisely identical. Therefore, we have also tested the robustness of our scheme to deviations 
of the individual systems from identicality. To this end, we replace F(xi) in Eq. (|30p by 



Fi(xi) 



-X i2 - X i3 

xu + 0.165a* 2 (l + A5 4 ) 
0.2 + (so. - I0)x a 



(30) 



where for each node i the parameter 6i is chosen randomly with uniform density in the interval [—1, 1]. We then repeat 
our original experiment. The parameter A characterizes the degree of non-identicality of the node dynamical systems. 
In Fig. [Til we show how the time averaged errors < E Xl > and < Ea > depend on A, for a network of N — 5 nodes, 
a = 0.2. As can be seen, the more non-identical the systems are, the worse our strategy performs. However, for small 
values of A, for example, for A < 0.03, the synchronization and identification errors are respectively less than 0.5% 
and 10%, i.e., < E Xl >< 5 x 10~ 3 and < Ea 0.1, thus indicating that good results may still be obtained when 
the coupled systems deviate from being precisely identical. 




FIG. 11: The time averaged errors < E X1 > and < Ea > versus the non-identicality parameter A, for a network of N = 5 
nodes, a = 0.2, Ti = 1, 6» = 0.6, v = 1/(ST S ), uj max = 2 x 10~ 4 , uj min = Lu max /2. 



D. Sensitivity to errors in the values of Aim 

The results presented previously rely on the assumption that the coupling with the maestro is exactly known at the 
network systems. Here we investigate a situation where the assumed known coupling with the maestro has a small 
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error, and we study the effects of this error on our strategy. To this aim we replace Eq. ([4]) by 

N 

ii(t) = F(xi(t))+r{r<-[— - 22 Ai j H(x i (t-T ij )) + (l-a)AiMH(xi(t-Ti M ))]} i i = 1, 2, N, , (31) 

where Aim = Mm^X + with 6\ is randomly chosen in [—1, 1], and the parameter S characterizes the degree to 
which the real coupling with the maestro Aim deviates from the coupling, AiM, assumed at node i. As can be seen 
from Fig. I12[ the larger the error, the worse our strategy performs. However, for small values of S, for example, for 
H < 0.03, the synchronization and identification errors are respectively less than 1% and 15%, i.e., < E Xl >< 10~ 2 
and < Ea >^ 0.15, thus indicating that good results may still be obtained provided that the coupling with the 
maestro is affected by a not too large error. 
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FIG. 12: The time averaged errors < E X1 > and < Ea > versus the mismatch parameter H, for a network of N = 5 nodes, 
a = 0.2, r x = 1, 6 = 0.6, v = 1/(8T S ), uj max = 2 x 10~ 4 , u mm = aw/2. 

We have also investigated the dependence of our strategy on the parameter a for the case in which small mismatches 
in the coupling with the maestro are present, i.e., S = 0.02. This is shown by the dotted line in Fig. [9] As can be 
seen, for such a case, our strategy works in a bounded interval of values of the parameter a, i.e., 0.1 < a < 0.45, 
indicating that under these conditions, the best setup is realized when synchronization depends on both the influence 
of the maestro, as well as that of the orchestra. 

IV. CONCLUSION 

In this paper we have dealt with the problem of identifying the time-evolutions of the couplings (assumed unknown) 
associated with the links of a network by making use of the phenomenon of chaos synchronization. We considered 
networks with time varying strengths of the couplings (due the effects of external unpredictable factors), as well as 
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different coupling delays affecting the communication between the nodes. 

We proposed an adaptive strategy based on the concept of a 'potential' that the network systems seek to minimize 
in order to identify the time evolutions of the network couplings. We considered fully connected networks of N nodes 
and showed how our strategy could be successfully applied to identify the N — 1 unknown couplings at each of the 
network nodes. We emphasize that, we achieved this goal by using very little information exchanged at the network 
nodes: namely the coupling was realized in only one state variable, and the only available information at each node at 
time t was an aggregate incoming signal from the other nodes, Ti(t) (compare for instance our approach with others 
previously reported in the literature [30I I31I]). Here we took advantage both of the phenomenon of synchronization 
of chaos and of the temporal signal diversity of the signal r-i(t) due to the presence of different time delays over the 
network links. 

We proposed two alternative versions of our strategy. In the first one, synchronization had to be achieved at each 
node (as well as the identification of the couplings over the network links) , based solely on the information received 
from the other nodes in the network. For this case, we showed that we were able to correctly identify the evolution 
of the couplings in the case of small networks, i.e., a network of N = 3 nodes, whose time evolution was sufficiently 
slow. 

In order to extend and improve our results, we introduced a second version of our strategy, where we made use 
of an additional node, termed the maestro, having the function of maintaining network synchronization. We made 
the assumption that the coupling with the maestro was known at the other network nodes and that the maestro was 
not affected by the dynamics of the other nodes in the network. We showed that such an arrangement could lead to 
much better results than in the case where the maestro was absent. In particular, we were able to identify the time 
evolution of networks that were larger and more rapidly evolving than was possible in the other case. Other possible 
setups of our strategy are also possible, e.g., a selective choice of the nodes to be connected to the maestro (pinning 
control) . 

V. APPENDIX 

In our experiment in Figs. 21 and [5] (Version 1, 6 = 0.3) we observed that we were able to synchronize the 
network on an approximate synchronous evolution obeying (J3]), but we were not able to correctly identify the time 
evolutions of the Ajj(t)'s. In what follows, we give an explanation for this finding. Our explanation is that, when 
the Ty are not large enough, it is possible to obtain approximate (not exact) synchronization for a range of values of 
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the Aij Aij. Although this synchronization is not exact, it is still good enough to imply that the minimum of the 
potential is so broad as to make our procedure ineffective. To illustrate this, assume the network systems evolve 
synchronously, i.e., x\{t) — X2(t) = xat(£) = x(t), and that the function H(x(t)) and x(t) are infinitely differentiable. 
Expanding H(x(t — Ty)) in a Mac-Laurin series in the parameter , we have 

H{ X {t-r ij )) = ±{- nj T^ d ^^. (32) 

n=0 

If this expansion of the function H(x(t — r)) has a finite radius of convergence, then (|32[) is valid under the assumption 
that Tjj is small enough. 

Our strategy seeks to minimize the potential at each node i, where VE^ may be thought of as a sliding average 
of the mean squared synchronization error at node i, and the synchronization error is 

N 

Yi(t) = n(t) -Y,^(t)H{ Xl {t - nj )). (33) 
Assuming Xi(t) = Xj(t) for all i and j and making use of (|32p . we obtain 

™=f: { ~T <flH lf )) £ [^w-^wiv (34) 

n=0 ' j'=lj'#i 

By construction, the Aij(t) change slowly with time (on the time scale ^ _1 ), as do the Aij(t) (on the time scale T n ). 
On the other hand, x(t), and hence - ^[^^ change rapidly with time. Thus to make the sliding average of Yi(t) 2 as 
small as possible, we need 

Aij(t) = A ij {t), (35) 

giving zero for the sliding average. However, if N is large and the Ty are not too big, we can achieve values of 
that are very small, even if (|3"5|) is far from being satisfied. To see this we note that if the delays ry are too small, 
for some sufficiently large n, say n* , the summation over n may be truncated with only a small resulting error in Yi. 
Assuming this to have been done, and setting each term in the sum over n in (|34[) equal to zero for each n < n* , we 
obtain n*N equations for the N(N — 1) unknowns A^. If the number of unknowns is greater than n*, then we can 
satisfy all the n* equations for a continuum of Aij values not satisfying (|35p . and the potential will then be near 
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zero. 
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